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Experimental moments of the nucleon structure function F2 
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Experimental data on the F% structure functions of the proton and deuteron, including recent results from 
CLAS at Jefferson Lab, have been used to construct their n < 12 moments. A comprehensive analysis of the 
moments in terms of the operator product expansion has been performed to separate the moments into leading 
and higher twist contributions. Particular attention was paid to the issue of nuclear corrections in the deuteron, 
when extracting the neutron moments from data. The difference between the proton and neutron moments was 
compared directly with lattice QCD simulations. Combining leading twist moments of the neutron and proton we 
found the d/u ratio at x — * 1 approaching 0, although the precision of the data did not allow to exclude the 1/5 
value. The higher twist components of the proton and neutron moments suggest that multi-parton correlations 
are isospin independent. 



The Operator Product Expansion (OPE) is a 
powerful tool in QCD which allows measurable 
moments of hadronic structure functions to be 
related to series expansions of the moments in 
terms of twists. The first term in the series, cor- 
responding to Leading Twist (LT), reflects the 
physics of asymptotic freedom, and is determined 
by single-parton distributions in the hadron. Sub- 
sequent terms in the series, or Higher Twists 
(HT), describe interactions between partons, or 
multi-parton correlations. The determination of 
the HTs is considerably more challenging, both 
experimentally and theoretically. 

We analyzed data on experimentally extracted 
moments of the proton and deuteron structure 
functions F2 JTJ to separate LT and HT terms. 
The n-th moment of the F2 structure function, 
including LT and HT contributions, can be writ- 
ten: 



M n (Q 2 ) = LT n (a s ) 
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where LT n is the leading, twist-2 moment, as is 
the running coupling constant, /i 2 is an arbitrary 
scale (taken to be 10 (GeV/c) 2 ), a T n is the matrix 
element of corresponding QCD operators, 7^ is 
the anomalous dimension and r is the order of 
the twist. 

The separation of the LT from the complete 
series is to some extent dependent on the or- 
der to which one calculates the LT Q 2 -evolution. 
In Fig. n we compare the n = 8 moment 
of the LT term calculated at fixed orders in 
pQCD: Leading Order (LO), Next-to- Leading Or- 
der (NLO) and Next-to-Next-to-Leading Order 
(NNLO); and by using resummation of soft gluon 
emission (SCR) 0: Leading Log (LL) and Next- 
to-Leading Log (NLL). The fixed order calcula- 
tions appear to converge at NLO 1 . However, at 
fixed pQCD order the logarithmic precision of the 
LT term deteriorates the closer one gets to x = 1 . 
Applying the SGR PJ we can improve the accu- 
racy of the LT term in the large-x region. The 
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iThe difference between NLO and NNLO is small with 
respect to uncertainties of the data 
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resummed LT calculated at NLL deviates signif- 
icantly from the LL result, similarly to the NLO 
and LO at fixed order. Unfortunately, next-to- 
next-to-leading log (NNLL) calculations are not 
yet available to confirm that also resummed mo- 
ments converge at second order. Nevertheless, the 
LT term is calculated to the best accuracy cur- 
rently available. 
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Figure 1. n = 8 LT moment calculated to dif- 
ferent pQCD accuracy: thin solid line - LO, thin 
dashed line - NLO, thin dotted line - NNLO, thick 
solid line - LL thick dashed line NLL. 



The extracted LT components of the proton 
and deuteron moments can be combined to form 
moments of the neutron Fi structure function. In 
the nuclear Impulse Approximation (IA), the nu- 
clear structure function can be written as a convo- 
lution of the nucleon structure function and a nu- 
cleon distribution function, f D , in the deuteron. 
In moment space this translates into a product 
of moments, so that the neutron moments can be 



obtained from: 

1M D ((1 2 ) 

M^m= n N { n ~*ff(Q 2 ). (2) 

where M?, M™ and are the proton, neutron 
and deuteron moments, respectively, and is 
the moment of the function f D . The distribution 
function f D was calculated from various deuteron 
wave functions 3 . 

The extracted LT proton and neutron moments 
can be combined to form Non-Singlet (NS) mo- 
ments of the nucleon F2 structure function, which 
can then be compared to lattice QCD simulations. 
A comparison of the extracted moments with re- 
cent lattice results from several groups is shown in 
Fig. |2 While a linear extrapolation of the lattice 
results to the physical pion mass overestimates 
our data significantly for n = 2 and 4, the re- 
sults extrapolated using chiral effective theory 
agree very well with our data. The data for higher 
moments are also of high precision, and it would 
be of considerable interest to compare these with 
higher lattice moments, especially since the ef- 
fects of chiral loops are expected to be suppressed 
in the large-n (large-x) domain. 

Another interesting result that can be obtained 
from the proton and neutron moments is related 
to the behavior of u and d quarks in the proton 
in the x — > 1 limit. At leading twist, the d/u 
ratio at large x can be extracted directly from the 
ratio of the neutron to proton structure functions 
F2 I F®, which is in turn related to the ratio of the 
moments M^/M? for large n. Indeed, for n 3> 1 
one finds that M£/M£(n -> 00) = F?/F%(x -> 
1). 

The n/p structure function moment ratios are 
shown in Fig. [3J Also indicated on the vertical 
axes are model predictions for the x — > 1 limits, 
namely the "standard" 1 /4 value used in most of 
parton distribution fits (corresponding to a van- 
ishing d/u ratio), and the value 3/7 expected from 
helicity conservation model (see Ref. 0). The 
trend of our data is towards the lower value of the 
model predictions with increasing n, although the 
precision of the data does not exclude the higher 
value. 

Indeed, at large n, both the data and the 
theoretical framework become more problematic, 
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Figure 2. Moments of the non-singlet F2 struc- 
ture function compared with several lattice QCD 
simulations: filled triangles - present analysis, 
open circles - lattice simulations from Rcf. 4 with 
linear extrapolation, open squares - lattice sim- 
ulations from Ref.jS] with linear extrapolation, 
open crosses - lattice simulations from Ref.jS] 
with chiral extrapolation from Ref.[5]. 



making it more difficult to distinguish between 
the different hypotheses. From Fig. one can 
also see the impact of the nuclear corrections on 
the deuteron moments. This is particularly evi- 
dent for large n: for n = 12, which corresponds 
to x values of around 0.75, the nuclear correction 
introduces a factor of ~ 2. 

Once the LT contribution to moments is deter- 
mined, one can then study the isospin dependence 
of the HT contribution. The HTs provide impor- 
tant information about multiparton correlations 
inside the nucleon. We assume that final state in- 
teractions and meson exchange currents (in par- 
ticular, the 1/Q 2 components) is negligible above 
Q 2 = 1 (GeV/c) 2 , so that the same nuclear cor- 



Figure 3. Ratio of neutron to proton moments 
as a function of n: filled triangles - this analy- 
sis, open circles - ratio without nuclear correc- 
tions applied. The solid (dashed) line indicates 
the scenario where d/u — ► 0(1/5) for x — ► 1 [7j. 



rections can be applied to the HTs as for the LT, 
Eq. [3 The total HT contributions to proton mo- 
ments and the corrected deuteron moments can 
then be compared, Fig. 01 This comparison indi- 
cates that the total HT contribution is indepen- 
dent of isospin. The isovector combinations p — n 
of structure functions F2 should therefore be free 
of HT contributions, within the presented uncer- 
tainties. 

In summary, we have analyzed experimental 
data on proton and deuteron F2 structure func- 
tion in order to extract their moments, and per- 
formed an OPE analysis to separate leading and 
higher twist contributions. By combining proton 
and deuteron moments and applying nuclear cor- 
rections, we extracted moments of the neutron F2 
structure function, paying particular attention to 
the issue of nuclear effects in the deuteron, which 
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Figure 4. Total higher twist contribution to 
the proton (filled triangles) and deuteron (open 
squares) moments at Q 2 = 2 (GeV/c) 2 . 



sociated with the pion cloud [H], but un- 
derestimate the lattice results when linearly 
extrapolations are used; 

• the total contribution of HTs is found to 
be isospin independent, which implies that 
in the isovector combination p — n of Fi 
structure functions the HTs are consistent 
with zero. 
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are increasingly important for higher moments. 

The LT part of the non-singlet moments were 
then obtained and related to lattice QCD mo- 
ments available for n = 2 and 4. For large n, the 
ratio of the neutron to proton moments can be 
related to the ratio of u and d quark contribu- 
tions in the proton in the x — > 1 limit. The HT 
contribution is related to the physics beyond the 
asymptotically free regime — namely, multipar- 
ton correlations. The results of our analysis can 
be summarized as follows: 

• the ratio of neutron to proton moments is 
consistent with F^ jF® — ► 1/4 as x — > 1 , al- 
though one cannot exclude the higher value 
of 3/7 suggested by helicity conservation ar- 
guments; 

• the non-singlet moments are in excellent 
agreement with the lattice data |4I5| . if 
these are extrapolated to physical quark 
masses taking into account chiral loops as- 



